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We consider the Standard Model with right-handed neutrinos to explain the masses of active 
neutrinos by the seesaw mechanism. Since active neutrinos as well as heavy neutral leptons are 
Majorana fermions in this case, the lepton number violating process can be induced. We discuss the 
inverse neutrinoless double beta decay e~e~ IV~W~ in the framework of the seesaw mechanism 
and its detectability at future colliders. It is shown that the cross section can be 17 fb for y/s = 3 
TeV even with the stringent constraint from the neutrinoless double beta decays if three (or more) 
right-handed neutrinos exist. In such a case, the future e“e“ colliders can test lepton number 
violation mediated by a right-handed neutrino lighter than about 10 TeV. 


I. INTRODUCTION 

The origin of neutrino masses is one of the most im¬ 
portant questions in particle physics at present. Various 
oscillation experiments have provided the mass squared 
differences and mixing angles of active neutrinos very pre¬ 
cisely [I|; however, the mass ordering, violation of CP 
symmetry and fundamental property of massive neutri¬ 
nos (i.e., Majorana or Dirac particles) are still unknown. 
The simplest way to explain the neutrino masses is to 
add right-handed neutrinos to the Standard Model (SM). 
The smallness of active neutrino masses can be explained 
by the large Majorana masses of right -handed neutrinos 
thanks to the seesaw mechanism [2 t81- 

In the seesaw mechanism, the mass eigenstates of neu¬ 
trinos are three active neutrinos which have tiny masses 
observed in oscillation experiments, and heavy neutral 
leptons (HNLs) which are almost identical to the right- 
handed states. The HNL masses are determined by 
Majorana masses and independent from the electroweak 
Higgs mechanism. They can take arbitrary values as long 
as they are sufficiently heavy to realize the seesaw mecha¬ 
nism. This is the reason why such HNLs can cause inter¬ 
esting phenomena in various aspects of particle physics 
and cosmology. 

Right-handed neutrinos can explain the baryon asym¬ 
metry of the Universe. The well-known scenario is the 
canonical leptogenesis [9| in which they need to be heav¬ 
ier than 0(10®) GeV [l^ (or 0(10®) GeV when the non- 
thermal production is realized [il|). The resonant lepto- 
genesis with quasidegenerate right-handed neutrinos [l^ 
can be effective with much smaller masses. Moreover, if 
we use the flavor oscillation of right-handed neutrinos, 
the required mass can be as small as 0(1) MeV p^l - [l5l| . 

In addition, a right-handed neutrino can play an im¬ 
portant role in astrophysics. It can be the dark matter 
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candidate with 0(1) keV mass M- This particle can also 
explain other phenomena, such as the pulsar kick [ia(for 
a review, see Ref. M)- Further, right-handed neutrinos 
with 0(0.Id GeV may be important for the supernova 
explosion [l9l |. 

Right-handed neutrinos are required by grand unified 
theories based on SO(IO) [1^ [2l| or larger symmetry 
groups. These particles can also accomplish the bottom- 
tau Yukawa unification [2^ which is necessary for real¬ 
istic grand unified theories. It has been discussed that 
they are favored to be lighter than O(IO^) GeV due to 
the argument of the electroweak naturalness [1^ . 

In this way, right-handed neutrinos in the seesaw mech¬ 
anism are well-motivated particles beyond the SM. The 
direct test of HNLs is possible if they are light and 
their mixings are large. So far various experiments have 
been performed to search for HNLs (for example, see re¬ 
views [23 - 1^ ). 

The noble consequence of the seesaw mechanism is that 
active neutrinos and HNLs are Majorana particles, and 
hence the lepton number violating processes can appear. 
A well-known example is the neutrinoless double beta 
(0^/3/3) decay {Z, A) ^ {Z -\- 2, A) -j- 2e“, which violates 
the lepton number by two units [28|. 

In this paper, we consider another example of such 
processes, 

e“e"^IU“lU“, (1) 

which is called the “inverse neutrinoless double beta de¬ 
cay” (see Fig. m, which cannot also happen in the 
SM. The e“e“ collision at high energy is one attractive 
option of the future lepton colliders such as the Interna¬ 
tional Linear Gollider (ILC [s^) and the Compact Linear 
Collider (CLIC [3l|). Various aspects of the process dD) 
have been studied so far [H, Isl - ldOl . 

The advantages of the process (HD over the OvfJfJ decay 
for the test of the lepton number violation are (i) the sig¬ 
nal event is clean and (ii) the prediction is free from the 
uncertainty in the nuclear matrix element. Furthermore, 
we should mention that these processes are complemen¬ 
tary tests for the lepton number violation. As we will 
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FIG. 1: Feynman diagrams for e e IV IV 


show, HNLs can give a destructive contribution to the 
decay so that no signal could be observed. Even in 
this case, the inverse Qvj3j3 decay could be observed since 
the energy scales of these processes are different. 

In this paper, we shall revisit the inverse Oz//3/3 decay 
paying special attention to the following points: (i) We 
study the process in the concrete model, i.e., the SM 
extended by A/" right-handed neutrinos with the seesaw 
mechanism, and try to derive the theoretical prediction 
being specific to the model, (ii) We take fully into ac¬ 
count the interference effects between active neutrinos 
and HNLs for both the Qvfili and the inverse Oz//3/3 pro¬ 
cesses. (iii) We estimate the upper bounds on the cross 
section depending on the number of right-handed neu¬ 
trinos. We use corrected cross section (for the difference 
from previous studies, see the Appendix). In addition 
we take into account the possibilities of the fine-tuning 
in parameters of the model, which have not been studied 
thoroughly. 

We find that the maximal cross section can be 0.47 
fb (17 fb) for the center-of-mass energy ^/s = 500 GeV 
(3 TeV) avoiding the stringent constraint from the 
decay if there are three right-handed neutrinos. This is 
a contrast to the previous results. Therefore, the process 
0 can be a good target of the future lepton colliders. 

The rest of this paper is organized as follows. In Sec. 
im we define the model parameters and list existing con¬ 
straints on them. In Sec Ell we derive upper bounds on 
the cross sections of the process © in the cases with one, 
two, three right-handed neutrinos in turn. Finally, Sec. 
m is devoted to conclusions. We add an Appendix to 
present the definitions of variables in the cross section of 
e-e- ^ W-W-. 


II. STANDARD MODEL WITH 
RIGHT-HANDED NEUTRINOS 

Let us first explain the framework of the present anal¬ 
ysis. We consider the SM extended by right-handed neu¬ 
trinos lyji with Lagrangian 

= >CsM + 

— Mj_ 

— FalLa^VRI - + h.C. . ( 2 ) 

Here we shall assume = lAa/K'h) <C Mj {a = 

e, /r, r; / = 1, • • • ,J\f) in order to realize the seesaw mech¬ 


anism. In this case the mass eigenstates of neutrinos are 
three active neutrinos Vi (i = 1,2,3) with masses rrii and 
Af HNLs Ni with masses ~ Mj. The mass ordering of 
HNLs can be chosen as Mi < M 2 < • • • < Mj^ without 
loss of generality. The mixing of neutrinos in the charged 
current interactions is then written as 

3 At 

ULol = ^ Uo^iUi + ^ QoliNj , (3) 

7^1 

where U is the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix of active neutrinos, while Qai 
represent the mixings of HNLs obtained by the diagonal- 
ization of the (3 -I- W) x (3 -I- Af) neutrino mass matrix 
including radiative corrections. 

In this model, the masses and mixings of active neu¬ 
trinos and HNLs must satisfy the following relations: 

3 At 

0 = ^ ^ e^iQ^iMi. (4) 

i 1=1 

These relations hold if the Majorana mass terms for the 

left-handed neutrinos are absent at the tree level in the 
canonical seesaw mechanism discussed here. We should 
note that such Majorana masses are induced by radiative 
corrections [4l|, and would alter the relation (g]) des¬ 
perately. We assume that such corrections are sufficiently 
suppressed by artificial fine-tuning. 

For the inverse QvjSp decay (H]), only the component 
with a = P = e. 


O = m:ff + ^02,Mz, (5) 

I 

is relevant. We call this equation “the seesaw relation” 
from now on. We dehne by 

3 

Ks = "mr ■ ( 6 ) 

i=l 

In the considered model both active neutrinos and 
HNLs are Majorana fermions, and then they induce 0^/3(3 
decay. The half-life of the decay is expressed as 

= (7) 

Here and hereafter, we use the notation and the results 
for the 0;//3/3 decay given in Ref. [i^. The effective neu¬ 
trino mass in the considered model is given by 

At 

meff = ^ QIj Mi fi, (8) 

7=1 

Here the first term represents the contribution from the 
active neutrinos which is given by Eq. ([6]), and the 
second term denotes the contributions from HNLs, in 
which the suppression of the nuclear matrix element for 
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M/ ^ 200 MeV is taken into account by the 

function // 


fi 


{P^) 

(p2) + Ml ■ 


( 9 ) 


As for the {)vl3j3 decay of ^^®Xe, the lower bound on the 
half-life is T 1/2 > 3.4 x 10^® yr at 90% C.L. [d^, which 
can be translated into the upper bound on the effective 
neutrino mass as 


ImeffI < = (0.185- 0.276) eV, (10) 


where we have taken into account the uncertainties esti¬ 
mated in [ 4 ^. On the other hand, the half-life bound of 
^®Ge is ri /2 > 3.0 X lO^® yr at 90% C.L. [H. In this 
case the bound on |meff| is 

ImeffI < = (0.213 - 0.308) eV. (11) 

We then apply the bound 113 throughout this analysis 
in order to make the conservative analysis. 

Furthermore, the mix ing s of HNLs are constrained by 
various experiments [2#4^ . In the following analysis, we 
consider the case when the mass of the lightest HNL is 
^ 3 GeV in order to avoid the stringent constraints 
from the search experiments by using the decays of tt, 
K and D mesons. Even in this case, there is an upper 
limit on |0e/P in the region Mj < mz from the search 
for HNLs at the Large Electron-Positron collider (LEP) 
experiment by the decay of the Z boson . In the mass 
range Mj ~ 6-50 GeV, the bound is given as |0e/P < 
2.1 X 10“® at 95% C.L. There is also an upper bound 
on the mixing angle which comes from the electroweak 
precision tests. The recent bound at 90 % C.L. [d^ is 


|0ep = ^ |0e/P < |0e||w = 2-1 X lO'^ . (12) 

I 

These are the upper bounds on the mixings which are 
important for the analysis below. ^ 


III. INVERSE NEUTRINOLESS DOUBLE BETA 
DECAY e"e" W W 


Now, we are in the position to discuss the inverse neu¬ 
trinoless double beta decay O- The cross section of the 
process in the model with Af right-handed neutrinos is 
given by 


dajif 

dcos 0 


GlPw 

32tt 


|Atp Bt -|- |Aup Bu 


+ {-At A* -1- h.c.) Btu 


(13) 


Here the definitions of variables in this equation are pre¬ 
sented in the Appendix. Notice that only At and A„ 
depend on the parameters of HNLs, namely, Mi and 0e/. 

In the following analysis, we consider the detectability 
of e~e~ —>■ W~W~ by taking the center-of-mass energy 
as ^/s = 0.5, 1 and 3 TeV. We assume the integrated 
luminosity of 100 fb, which can be achieved by 0(10^) s 
run of the ILC or CLIC [d^ (the luminosities of the e“e“ 
colliders are expected to be the same order as e+e“ col¬ 
liders [sij). 

Before discussing the realistic cases, we shall consider 
the two extreme cases. The first one is the limit when 
all HNLs are sufficiently heavy (i.e., Mj ^ -^ 5 ) and only 
three active neutrinos effectively take part in the process. 
In this case the cross section is given by 


ao = 


3G%1 


‘eff I 


dTT 


= 0.96 X 10 


-18 


fb 



The cross section is determined solely by the effective 
neutrino mass rUgg of active neutrinos in the Oi^/S/S decay. 
It is seen that (Tq is too small to be accessible by future 
colliders. 

The other limit is that all HNLs are sufficiently light 
(i.e., Ml mw ^ y^). In this case it is found from the 
seesaw relation that 


1 


M 


At — - TOgg + y 0 gj Ml — 0 


7=1 


(15) 


and then A^ = 0, which results in the vanishing cross 
section. In other words, the seesaw relation m ensures 
the unitarity of the process o at high energies y/s 

Ml [ 13 . 

From now on, we will study the cross section of 
e~e~ —>■ W~W~ in the framework of the seesaw model 
with Af = I, 2, and 3 right-handed neutrino(s). Espe¬ 
cially, we shall discuss the impacts of various constraints 
on the masses and mixings of HNLs discussed in the pre¬ 
vious section. 


A. Case with one right-handed neutrino 

Let us first consider the SM with only one right-handed 
neutrino. In this case, there is one HNL Ni with mass 
Ml and mixing 0 q,i in addition to active neutrinos. This 
model is incompatible with the oscillation data since 
there appear two massless states among active neutri¬ 
nos. We shall, however, discuss it in order to make our 
point clear. 

One might think that Ni with a relatively large mix¬ 
ing allowed by Eq. m may give rise to a significant 


^ Our analysis can also be applied to processes — >■ W~W~ 

and e~W~W~ as in Ref. 1^ . The latter process is, 
however, strongly constrained by the /i —>• ey experiment [4 ^l48l| . 

IE/©:7©m/I <10-5 (90%C.L.). 


^ The prefactor of ao in Eq. m is three times larger than the 
estimations in Refs. [33J39|| (see the discussion in the Appendix). 
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FIG. 2: Cross section of e~e~ W~W~ in the model with 
one right-handed neutrino [M = 1) for ^/s = 3 TeV (the 
red solid line), 1 TeV (the green dashed line), and 0.5 TeV 
(the blue dotted-dashed line). The horizontal line (the black 
dotted line) represents no in (1141) . Here we take \mgg\ = 0.276 
eV. 


contribution to e“e“ —>■ W~W~. This is, however, not 
true in the considered case. The seesaw relation (0 tells 
that the mixing cannot be large but Oh = 
and then At is suppressed as 


This shows that the cross section becomes independent 
on the mixing of but is determined from the effective 
mass TOgjj and Mi. In Fig. [5] we show the cross section 
e-e- ^ W-W- with \m’^g\ = 0.276 eV. 

It is seen that the cross section approaches to uq IT^ 
for Ml ^ y/s due to the decoupling of Ni as mentioned 
before. On the other hand, when Mi <C mw-, the cross 
section behaves as 


CTl 


1 

Gtt 


(17) 


and it vanishes for Mi —>■ 0 as expected. 

In the model with Af = 1, therefore, the cross section 
is too small to be observed at future colliders. This is 
the direct consequence of the seesaw mechanism, namely, 
the mixing 0 ei is determined by from the seesaw 
relation 


(be., mi = 0 for the normal hierarchy and m 3 = 0 for 
the inverted hierarchy). 

In this case, by using the seesaw relation ([5|) and |0ep 
(na), the mixing angles of and N 2 are given by^ 


0 

0 


2 

el 


2 

e2 


M 2 |2 m: 


Ml + M 2 


|0e|^- 


eff 


Ml + M 2 


Mi |e^p_ "IC, 


MI + M 2 


MI + M 2 


(18) 

(19) 


We find that, when Mi <C M 2 and jm^fjl/Mi <C |0eP, 
the large mixing of Ni, Qh ~ | 0 e|EWi realized. 

This is because of the cancellation between the contribu¬ 
tions of and N 2 in the seesaw relation © [Hil. 
Some mechanism is desirable to stabilize such a fine- 
tuning, for example, by a discrete flavor symmetry. This 
issue is, however, beyond the scope of our a naly sis (sim¬ 
ilar situations have been discussed in Refs. [50,[51|). 

The function At is now 


At 


{M2 — Ml) Ml M2 I „ 1 2 

{t-Ml){t-Ml) 

[(^2 - Ml M 2 + Ml) t - Ml M|] 
t{t - Ml){t - Ml) 


( 20 ) 


It is then seen that, by neglecting the second term sup¬ 
pressed by m’^g, the hierarchical mass pattern Mi <C M 2 
enhances the matrix element of the process. In this limit. 
At is dominated by the contribution of iVi as 


At 


Ml |0eP 

t-Ml ’ 


( 21 ) 


which leads to the cross section 
G|s|0e|4 


02 ^ 


Stt 


H{ri ), 


( 22 ) 


where ri = Ml/s and 77(ri) is 77 ~ 6 ri for Mi <C mwi 
H ~ 2ri for mw “C Mi <C y/s, and H ~ l/(2ri) for 
Ml ^ y/s. By taking mw/y/s —t 0, it is approximately 
given by 


H{r) 


r(2 -I- 3r) 


1 

1 -I- r 


2r 


l + 2r 


■ In 



• (23) 


It is then found that H{ri) takes its maximal value 

i?(ri)|MAX = 0.206 (24) 

at ri = 0.542 (be.. Mi = 0.736i/s). If we applied the 
upper bound on |0eP given in Eq. (TT^ . the maximal 
cross section 0-2 = 17 fb could be obtained for yfs = 3 


B. Case with two right-handed neutrinos 

Next, we discuss the case with two right-handed neu¬ 
trinos, where there appear two HNLs and 7 V 2 . Notice 
that the lightest active neutrino in this case is massless 


® can be taken to be real and positive without loss of gener¬ 

ality. Here we consider the case when 0^2 is real and negative. 
The importance of the relative phases between the mixings has 
been discussed in Refs. l36l. l37ll . 
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FIG. 3: Cross section of e~e~ —>■ W~W~ in the model with 
two right-handed neutrinos [N = 2) for yfs = 3 TeV. We take 
M 2 = l.OlMi (the blue dotted-dashed line), M 2 = I.IM 2 (the 
green dashed line) and M 2 = IOOM 2 (the red solid line). Here 
we consider the inverted hierarchy case with = 0.276 eV 
and = 4.79 x 10“^ eV. 


TeV. Unfortunately, such a large cross section cannot be 
realized for the case with A/" = 2, because we have to take 
into account the severe constraint from the OvPP decay. 
The effective neutrino mass in this case is estimated as 


rUeS = 


ifi 


/ 2 ) 


Ml M 2 
Ml +M 2 


|0eP 



Ml + M 2 ) 


(25) 


The upper bound on |0ep is then given by Eq. (IT^ for 
Ml > M*, while the more stringent bound from the QvjSfi 
decay^ 


,0 ,2 , + |TO^ff|)MiM 2 

^ (p2) (M 2 - Ml) 


(26) 


is imposed for Mi < M* where M* is determined from 
the ratio Mi/M 2 as 


M* 


IQellw(p^) 

V M2) -b Im/ffl ■ 


(27) 


Here we have assumed that Mf^ M| ^ (p^)- 

By using the upper bound on the mixing angle listed 
above, we can estimate the maximal value of the cross 


^ When the masses of A^i and N 2 are sufficiently degenerate, we 
can avoid the constraint from the Oupp decay. In this case, how¬ 
ever, At in Eq. approaches to ca, which leads to the sup¬ 
pression of the cross section (see also Fig. EJ. 


FIG. 4: Cross section of e“e“ —^ W~W~ in the model with 
two right-handed neutrinos with Mi <C M2 for ^/s = 3 TeV 
(the red solid lines), 1 TeV (the green dashed lines), and 0.5 
TeV (the blue dotted-dashed lines). In each set of lines, the 
upper one is for the inverted hierarchy case with \m'/g\ = 
4.79 X 10~® eV and the lower one is for the normal hierarchy 
case with Im/gl = 3.66 x 10“® eV. We use = 0.276 and 
0.185 eV for the thick and thin lines. 


section of e~e~ —>■ W~W~. In Fig. [3]we show how this 
maximal value depends on the mass ratio when y/s = 3 
TeV. Here we choose M 2 /M 1 = 1.01, 1.1 and 100. It can 
be seen that the largest cross section is obtained by taking 
M 2 /Ml —>■ 00 . (Note that the difference in the cross 
section between M 2 /Mi = 00 and 100 is very small.) 

We then show in Fig. S] the upper bound on the cross 
section by taking y/s = 0.5, 1 and 3 TeV. Notice that, 
when Ml < M* = 0(10^) GeV, the maximal cross sec¬ 
tion is given by 


(28) 


Stt 




The function ri Hiji) becomes constant for ri ^ 1. For 
Ml > M* the cross section falls as r//^. The maximal 
cross section is then given by 


0-2 < 8.9 x lO-^fb 


r^eS+lKsl y 

(0.276 -b 0.048) eV J 


f 0.178 GevV ( ^ \ 

v^/PT; V3^; 


(29) 


The maximal value changes depending on the mass 
hierarchy of active neutrinos. The cross section in the 
inverted hierarchy can be larger than that in the normal 
hierarchy since |mg(f| can be larger: < 3.66 x 10“^ 

eV for the normal hierarchy and Im/ffl < 4.79 x 10“^ eV 
for the inverted hierarchy. Here we have estimated these 
values by using the central values of the mass squared 
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differences and the mixing angles of active neutrinos in 
Ref. [I| and by varying the possible values of the CP 
violating phases in the PMNS matrix. Notice that, as 
shown in Fig. 21 there exists an uncertainty in the up¬ 
per bound on the cross section from the nuclear matrix 
element shown in Eq. (uni). 

The upper bound on the cross section scales as and 
hence the collisions of electrons with higher energies are 
desired. In fact, the possible way to observe the inverse 
OvPP decay is to realize ^ = 3 TeV at CLIC with the 
integrated luminosity larger than 0(10^) fb. In such a 
case, the number of the signal event can be larger than 
unity. 


C. Case with three right-handed neutrinos 

Finally, let us discuss the case with three right-handed 
neutrinos {Af = 3). We will show that this case allows 
the larger cross section than the previous cases. Our 
basic idea is the following: Let consider three HNLs with 
hierarchical masses Mi <C M 2 <C M 3 . Then, the large 
cross section of e~e~ —>■ W~W~ can be induced by N 2 
with mass M 2 ^ -v/s and mixing | 0 e 2 p — lOellw ©■ 
The seesaw relation can be realized even with the large 
mixing 1 0 e 2 P because of the cancellation between N 2 and 
N 3 , which is similar to the case with Af = 2. In addition, 
the constraint from the QvfAji decay can be avoided by 
the cancellation in rUeS between N 2 and A^i. This is the 
reason why A/" > 3 is required for having the large cross 
section of the inverse QvPP decay in the framework of the 
seesaw mechanism. 

Now, we discuss this point in detail. By using the 
seesaw relation and the effective mass (|S]), we may 
express the mixing angles of and N 3 in terms of that 
of N 2 '. 


0 


2 

el 


-I 

(/i - h) Ml 


(/2 


/a) M 2 0^2 



FIG. 5: Cross section of e“e“ —^ W~IV~ in the model with 
three right-handed neutrinos (A/ = 3) with Mi <C M 2 <C M 3 
for = 3 TeV (the red solid line), 1 TeV (the green dashed 
line), and 0.5 TeV (the blue dotted-dashed line). 


the masses of three HNLs are hierarchical, namely Mi <C 
M2 M3. In this case the mixing angles become 


0'i 


0e3 ^ 


_^02 


M 2 


e2 5 


- — 02 
M 3 


(33) 

(34) 


and hence |0elp, |0e3p < | 0 e 2 p, so | 0 e 2 p |0ep, 

which results in 


M2|0e|2 

t-Mi ’ 

In this limit the cross section is given by 


(35) 


- TOeff + (1 - /a) m^s : (30) 


-I 


-b TOeff - (I - /l) TOgfi , 


(31) 


It is then found that At is 

^ M 2 (M| - Ml) {Ml - M|) {t + (p2)) 02^ 
* (M| + (p2)) [t - Ml) {t - Ml) {t - Mi) 

TOeff (Ml + (p2)) (M| -b (p2)) 

(p2) (^t-Ml){t-Mi) 

I Ke^lMi {t + (p^)) 
^{p^)t{t-Ml)it-Mi)- ^ 


0-3 ~ 


Gls\e. 

Stt 


ft |4 


-H{r 2 ), 


(36) 


where r 2 = Ml/s. The result is similar to Eq. (l22ll for 
N =2. 

Thus, the upper bound on the cross section can be 
obtained by the maximal value of the mixing angle, which 
is shown in Fig. [5] Here we take Mi = 3 GeV to avoid 
the stringent constraint on |0eiP from the direct search 
experiments. Notice that the bound on |0e|2 comes from 
the search at LEP for M 2 ^ rnz, and comes from 
the electroweak precision test m for M 2 P rnz. Since 
the former bound is so stringent, the cross section is too 
small to be observed for that mass range. The maximal 
cross section is obtained when M 2 — 0.736-v/i as 


By neglecting the second and third terms suppressed by 
TOeff and TOgjj, the matrix element can be maximal when 


(73 < 17 fb 


IQeP 

3TeV; \\eegw 


2 


(37) 
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FIG. 6: Sensitivity limit of |0e2p from e“e“ —>■ W~W~ in 
the model with three right-handed neutrinos (Af = 3) for 
y/s — 3 TeV (the red solid line), 1 TeV (the green dashed 
line), and 0.5 TeV (the blue dotted-dashed line) . The black 
dotted lines show the upper bounds from Refs. [4^ and 
The black long-dashed line shows the upper bound from the 
perturbativity of Yukawa couplings when M 3 = 10® GeV. 

Compared with Eq. (HH), the maximal value of the cross 
section in the Af = 3 case can be much larger than the 
previous cases. 

In Fig. ini we show the sensitivity limit of |0e2p to¬ 
gether with the current upper limits (for the perturba¬ 
tivity bound |0e2p < /{M 2 M 3 ), see Ref. [^). 

Interestingly, the inverse Ou/S/S decay process can probe 
HNL {N 2 in this case) of which the mass is much heavier 
than the center of mass energy (M 2 » i/i), because it 
is induced as the virtual effect of HNLs. This is one of 
the most important advantages of e~e~ —?> W~W~ in 
searching HNLs [H, . 

It should be noted that the process (HD is induced by 
the collision of two electrons with negative helicity. Thus, 
this reaction will be turned off simply using one of the 
electron beams with positive helicity |32l| . This can offer 
the important test for this process by using the polarized 
beam which is specific to the linear colliders. 

We should note that the large cross section shown in 
Eq. ()37p is possible when the HNL with Mi <C M 2 
and 0gi ~ —{M 1 /M 2 ) QI 2 exists in order to avoid the 
stringent constraint from the Qvfip decay. If this is the 
case, such a light HNL Ni is also a good target for the fu¬ 
ture experiments. The searches b y u sing charmed meson 
decay e+e“ —?> uNi [ioL issjl and Z —^ vNi [ 5 ^ 

may provide the complementary test for the model dis¬ 
cussed in this analysis. 

Here, we have derived the sensitivity limit of the mix¬ 
ing by just counting the numbers of the signal event as¬ 
suming no background event. The pos sible background 
processes have been studied in , which may 


reduce the sensitivity. However, the usage of the polar¬ 
ized beam can make the cross section four times larger 
and the longer duration of the experiment can increase 
the number of events. These issues must be studied in 
detail to have the precise sensitivity limit; however, they 
are beyond the scope of this analysis. 

So far, we have discussed the case with the hierarchical 
mass pattern Mi <C M 2 M 3 with M 2 ^ y/s. When 
M 2 ^ yfs, N 2 decouples from the process and the cross 
section becomes small. Even in this case, when Mi ~ 
yfs <C M 2 M 3 , the situation becomes the same as the 
case with Af = 2. Then, the maximal cross section is 
given by Eq. (I^ . 

Note that Eq. ((36D is independent from meg. This 
means that the inverse 0^'/3/3 decay can happen even if 
the the Ov^fi decay could not be observed. This is be¬ 
cause the energy scales of the former process (i/s ^ 1 
TeV) and the latter one {yj (p^) ^0.1 GeV) are much 
different. Therefore, the 0 i 2 / 3/3 decay and the inverse one 
are regarded as complementary tests for the violation of 
lepton number in nature. 

Finally, we comment on the case with Af > 3. The 
maximal cross section is the same as that with Af = 3. 
This is because the severe constraints from the seesaw 
relation and the Qvj3fi decay can be avoided as described 
above, but the inevitable upper bound (IT^ on | 0 eP, 
which is just the sum of | 0 e/p, gives the cross section 
as shown in Eq. (|37|. 


IV. CONCLUSIONS 

We have discussed the inverse Qvj3fi decay process 
e“e“ —>■ W~W~ in the SM with Af right-handed neu¬ 
trinos. We have estimated for the first time the max¬ 
imal cross sections of this process depending on Af by 
taking into account the seesaw relation on the mixings 
of active neutrinos and HNLs as well as the severe con¬ 
straints from direct search experiments, electroweak pre¬ 
cision measurements and 0^/3/3 decay. The maximal cross 
sections are realized in limited regions of the parameter 
space, and the cross sections are orders of magnitude 
smaller if there is no fine-tuning. 

For the Af = 1 case, the seesaw relation restricts the 
cross section to be smaller than 0 ( 10 “^®) fb, so it cannot 
be seen at future experiments. For the Af = 2 case, the 
cross section can be a 2 ^ 10“^ fb for y/s = 3 TeV. It 
may then be observable at CLIC, but it is impossible at 
ILC. For the Af = 3 case, the larger cross section as 0-3 = 
0.47 (17) fb for y/s = 0.5 (3) TeV can be obtained. In 
this case both ILC and CLIC can search the inverse 
decay, which would reveal the fate of the lepton number 
and the origin of the neutrino masses. 
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Appendix A: DEFINITION OF VARIABLES 


The range of |t| is given by (1 — 2 +/3w)s/2 > \t\ > 
{1 — 2rw — I3w)sl2 and then |t| > [rw + 0(r^)]. 
Therefore, mf/ltl <C 1 for the realistic range of ^/s. In 
this case. At is 


At 


™eff i QeJ 

t ^ t-Mf ' 

1=1 I 


(A5) 


In this Appendix we define the variables used in the 
cross section of e~e~ —>■ W~W~ in Eq. m- Here, we 
work in the center-of-mass frame and 0 is a scattering 
angle of one of the outgoing W~s. Mandelstam variables 
are defined as usual, 

t =-|(l-/3tv cos 6») + m^, (Al) 

M =-|(1+/3wcos6»)+ m^, (A2) 

where Pw = y/T^^4rw with rw = w^/s. The variables 
At and A^ are given by 




Uei 

t-ml 


A eljMi 
^ t-Mf^ 


Au — 




t—¥u ' 


(A3) 

(A4) 


We have found a difference in the definition of At from 
Refs. [ 3 ^, 113 • Here we have corrected the errors in the 
denominators of At as 2 t —>■ t. 


On the other hand, the variables Bt, Bu and Btu are 
given by 


Bt 

= (1 - 4rv^)f^ - 4rw(l - 2rw) sf 



+ ^{l-rw)rw 

(A6) 

Bu 

Bt\t=^u ’ 

(A7) 

Btu 

= {1 - Arw)tuP . 

(AS) 


A typo in Ref. [s^ is corrected in Eq. (IA6I) . 

Since we are considering y/s > 500 GeV, we can ap¬ 
proximate rriw “C s. We cannot, however, approximate 
|t|, since |f| ~ rn^js < at 0 ~ 0. Indeed, 

/ I 4 

(icos0^^^A = 2 -b C>(rw), (A9) 


so we have to take the limit rn^js —^ 0 after integra¬ 
tion when we consider the contribution of light neutrinos. 
This is the reason that the cross section in Eq. 
three times larger than that in previous studies 
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